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1. Introduction. 

It is well known that the principal chiral SU(N) x SU(N) model (PCM) in two Euclidean 
dimensions is similar to four dimensional pure SU(N) gauge theory in many respects Q. 

Recent numerical work provides evidence that Wilson loops in SU(N) gauge theory in 
two, three and four dimensions exhibit an infinite N phase transition as they are dilated 
from a small size to a large one; in the course of this dilation the eigenvalue distribution 
of the untraced Wilson loop unitary matrix expands from a small arc on the unit circle 
to encompassing the entire unit circle |2|, ||. There is further evidence that in the vicinity 
of the critical size and for eigenvalues close to -1, the eigenvalues behave in a universal 
manner controlled by two critical exponents of N, taking the values 1/2 and 3/4. The uni- 
versality class of this transition is believed to be that of a random multiplicative ensemble 
of unitary matrices. The transition was discovered by Durhuus and Olesen Q when they 
solved the Migdal-Makeenko || loop equations in two dimensional planar QCD. We refer 
to it therefore as the DO transition and to the universality class as the DO class. The 
multiplicative random matrix ensemble || can be axiomatized in the language of noncom- 
mutative probability |^]. In some sense it provides a generalization of the familiar law of 
large numbers, which is associated to the abelian case. The essence of the difference is not 
in that the law of large numbers is additive, while the new case is multiplicative; rather, 
the essential features making a difference are that one case is commutative and the other 
not. Various recent insights into the DO transition |8], |9[ |l0| indicate that an even deeper 
understanding of the transition might emerge. 
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In this work we present numerical results indicating that the PCM model has a similar 
transition, this time in the two point correlator matrix: when it is dilated it undergoes 
a DO transition. Thus, the analogy between the PCM and gauge theory is upheld and 
the universal character of the DO transition as a marker for the transition scale sepa- 
rating weakly from strongly interacting dynamics in theories based on group manifolds is 
expanded. 

2. Basic facts about the PCM. 

The basic degrees of freedom are g(x) G SU(N) where x E R 2 and the action is: 

N f 

S = - J d 2 xTrd^g{x)d^{x) (2.1) 

The large N limit is taken in the 't Hooft prescription, by keeping the coupling T fixed, 
which makes S of naive order N 2 . 

There is a global symmetry group SU{N)i x SU{N)r acting on g{x) by g{x) — > 
Lg(x)I$ where L,R € SU(N). If we eliminate one of the factors by a translation breaking 
"gauge choice" g(0) = 1, we are left with a global symmetry given by a single SU(N) 
acting on g(x) by conjugation and leaving g(0) = 1 unchanged. This single SU(N) is the 
"diagonal subgroup" of the global symmetry group. 

The model is asymptotically free and has a nonzero massgap. Using Bethe ansatz 
methods it was found that there are N — 1 particle states with masses given by 

sinf— ) 

M r = M — i£f, l<r<JV-l (2.2) 
sin(f ) 

Under the diagonal symmetry group, the states corresponding to the r-th mass are a 
multiplet transforming as an r component antisymmetric tensor. This has been verified 
numerically for N = 6 in [|ll]]. In consequence it is believed that for any N and 1 < r < N— 1 
we would have 

(Xr(^(0)^(x)))~C P ^e- w -W (2 .3) 

for \x\M » 1. Here, Xr takes the trace in the r-antisymmetric representation, normalized 
by Xr (l) = ( N r ). 

The PCM can be formulated in terms of the currents J^(x) = g^ (x)d ll g(x). Note that 
the object g(0)g'(x) is determined by the via a path ordered exponent. In particular, 
also the action can be expressed in terms of the J fl (x). In addition, the J/i(x) obey local 
constraints because if viewed as gauge fields they are pure gauge. 

3. The average characteristic polynomial and the DO transition. 

It was found that a convenient observable to use in order to locate the DO transition 
and identify the associated critical iV-exponents is provided by the average characteristic 
polynomial of the unitary matrix undergoing the transition. This matrix is the open Wilson 
loop in the gauge case while in the PCM model it is g(Q)g^ (x). 
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For gauge theory in two dimensions, for a non-self-intersecting loop, the unitary matrix 
can be written as the product of many such open Wilson loop matrices round small loops 
tessellating the spanning area of the loop. For a gauge theory in higher dimensions we 
choose an infinite two dimensional surface extending to infinity and containing the loop in 
question and repeat the argument for gauge fields tangent to the surface. The DO transition 
occurs because for large loops the correlations between well separated tessellating loops die 
out and the factors become effectively independent. In two dimensions any two distinct 
tessellating loops are uncorrelated. Note that the "number" of factors in the product is 
proportional to the area enclosed by the loop; this product is not to be identified with the 
linearly ordered product associated with the path ordered exponential. 

For the PCM, we split the segment connecting to x into small intervals, with J M (y) 
determining a multiplicative contribution for each element. Again, for large the corre- 
lations between well separated small intervals dies out, and a DO transition might occur. 
This time however, the "number" of factors is proportional to the linear distance \x\. Again, 
this random product is not to be identified with the product between the two matrices in 
the two-point correlation function. And, yet again, if we descend in dimensions (here to 
one dimension) any two distinct intervals are uncorrelated and the DO transition can be 
established by an exact analytical solution. 

In either case, gauge or PCM, if W is the unitary matrix given by the product of 
independent factors, the average characteristic polynomial is defined in the same manner: 

N 

Q N (z, A) = <det(z - W)) = Y,z N - r (-l) r ( X r(W)) (3.1) 

r=0 

Here A is a continuous parameter which has been kept finite by tuning the distribution 
of the individual factors in the product to almost always produce something finite after 
multiplying the many factors. In a sense, A would measure the area in the gauge case and 
the distance in the PCM case. In general, one can view A as a parameter quantifying the 
amount of dilation of the object relative to a fixed standard. 

In the simple case of two dimensional gauge theory one can take A to be literally the 
area, made dimensionless with the aid of the coupling constant, and the DO transition can 
be seen when one uses 

(Xr(W)) = 

For the PCM, we expect, 



(Xr(W)) ~ C r 

We observe that the exponents in the gauge and PCM case are similar, with M\x\ playing 
the role of A/2. More precisely, for r = 0(1) or N — r = 0(1) and large N they match and 
for r ~ N they do no differ that much. Thus, if the coefficients C r vary less than the other 
factors, the singularity structure at a critical value of A in the infinite N limit would be the 
same. We conclude that the PCM model could also exhibit the DO transition. It would be 



N\ y r(N-r) 

e 2jv 



N 



1 1 sin(-&) 



(3.2) 



(3.3) 
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nice to confirm this by exact analytical methods, in particular since the numerical work we 
shall present indicates that the DO transition does occur in the PCM. The commonality 
between the gauge and PCM case exponents is related to a topic known as "Casimir 
scaling" . 



4. Regularization and renormalization. 

Until now we ignored renormalization. In both the gauge and PCM case one needs to 
renormalize and there is no guarantee that after that it makes sense to view W as a 
fluctuating unitary matrix of unit determinant. Formally, we can define a renormalized 
version of the W matrix in each case which does provide for such a view. The gauge case 
was analyzed in other work, so from now on we restrict ourselves to the PCM case. 

The extra regularization we put in is of the operator. The action is regularized with 
appropriate counter terms in some conventional manner. The net (formal) result is that 
we assume that we have a probability distribution generating "bare" configurations go(x). 

Let r > be a parameter of dimension length squared which means that tM 2 is a 
finite unitless parameter. We extend R 2 (space time) by an R + coordinated by r and define 
a g(x,r) determined by go(x) by: 

5 t (S)T) M^r) = l [S 2 5(a . jT) _ 5 V(x,r)] - ±Tr[d 2 g(x,r) - d 2 g\x,r)] (4.1) 

with the initial condition 

g(x,0)=g (x) (4.2) 

The evolution in r ensures that g(x,r) G SU(N) if g (x) G SU(N). For g{x,r) = e iB(x ^ 
with TrB(x,r) = and ||5|| << 1 the equation becomes 

23*2 »«»*(,, r) (4.3) 

This shows that for r > high momentum modes will get suppressed. g(x,r) depends 
nonlocally on go(y)', the amount of nonlocality is characterized by the unitless finite quantity 
tM 2 , which is parametrically small relative to one, but nonzero. As a result, the correlation 
function (g(x,r)g' (y,r)) will stay finite and nontrivial in the continuum limit so long as 
r > without needing the extraction of wave function renormalization factors. 

Hence, the SU(N) matrix W we shall compute the average characteristic polynomial 
of will depend on an extra parameter r > 0. The critical scale \x — y\ c will have some weak 
dependence on r. The point is that so long as tM 2 is parametrically small relative to unity, 
there will be some finite \x — y\ c where the transition occurs, but the exact value of that 
\x — y\ c is renormalization dependent. The relevant fact is that the DO transition occurs 
somewhere and there it marks a sharp boundary separating strong and weak couplings as 
defined from this particular two point function. 
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5. Lattice formulation. 



The continuum formulation ended up with three continuum coordinates. We now discretize 
each. The two dimensional space-time becomes a square lattice with unit lattice spacing, 
whose sites are again denoted by x, with x G Z 2 . ±/i denote unit vectors in the posi- 
tive/negative [i = l,2 direction. The r direction becomes a unit spaced line labeled by a 
coordinate t = 0,1,2,.... The variables that enter the action are g(x) € SU(N) and the 
r-dependence of the smeared variables is denoted by gt(x), with go(x) = g(x). The system 
is made finite by taking the x-space to be a torus lattice of equal side lengths, L. The ex- 
tent in the r direction is n. Only the space-time direction is stored in the computer, while 
the t dependence is iteratively computed. The g(x) obey periodic boundary conditions and 
gt(x) for £=1,2, n is determined by g(x) in a manner preserving this. 
The lattice action is: 

S L = -2Nb^ j m , r[g(x)g ] {x + [jl)\ (5.1) 

Averages with respect to 

Afl[[dg(x)}e-^ (5.2) 

X 

are denoted by (...) where M is a normalization constant ensuring (1) = 1. Thus, the two 
point function is 

G{x) = ±(Tr\g(0)J (*)]> (5-3) 

When N is taken to oo with b fixed (this is the 't Hooft large N limit) G{x) remains 
a finite nontrivial function of x. In particular, £q, defined below, measures the correlation 
in lattice units. 

2 _1 E^ 2 G(X) 

^ - 4 E x G(x) (5 ' 4) 
£g is independent of the wave function renormalization constant. It sets the scale for the 
system in units of the lattice spacing. The continuum limit is obtained when b is taken to 
infinity, which causes £g to diverge. 

It has been found that |llj] (a number further improved in subsequent work) that in 
the continuum one has 

M^ G = 0.991(1) (5.5) 

Here, M is the lattice equivalent of the continuum mass M = Mi from Equation fl2,2|) . We 
shall use £g as the scale setting quantity. 

The extrapolation to infinite b is greatly helped by using E, the link energy, as an 
expansion parameter: 

E = l-^(Tr[g(0)gHi)]) (5.6) 

As b — > oo, E — > 0. This scheme is equivalent to what is also known as "mean field 
improvement". In the range we shall be working and setting N to infinity, the following 
formula translates E into b to high enough accuracy: 

1 1 0.000545 0.00095 0.00043 
86 + 2W + p W~ + 6 5 ( } 
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The following formula is found to work within a few percent in the range we shall be 
working in, which is 11 < £g < 20. Again, we use the N = oo expression, but know that 



■ffz corrections are small. 



Cg = 0.991 



2-TT 

e 4 
16vr 



VE exp (- j (5.8) 



Using this formula as setting the scale for all dimensionful physical quantities, the approach 
to continuum can be tested to a level below one percent. 

The smearing operation is discretized as follows: We start from a configuration g(x) = 
go(x) and evolve with a fixed smearing parameter < / << 1 to a configuration gt(x). 
One smearing step takes us from gt{x) to gt+i{x). We first define Z t+ i(x) by: 

Z t+1 (x) = Y^lsitfrM* + - !] ( 5 - 9 ) 

Next we construct antihermitian traceless SU(N) matrices At+\(x) 

A t+1 (x) = Z t+1 (x) - z\ +l (x) - iTr(Z t+1 (x) - Z t f +1 (x)) = -A t f +1 (x) (5.10) 

and set 

L t+l {x) =exp[fA t +i(x)] (5.11) 
Finally, gt+i (x) is defined in terms of L t+ \ (x) by: 

9t+i(x) = g t (x)L t+1 (x) (5.12) 

This procedure is iterated until t = n is reached. Each iteration is performed for all x 
and only then one goes to the next iteration. 

For a given £q (or b) the parameter r is fixed keeping the quantity t/£ g unchanged; 
here, r = nf is in units of lattice spacing. Again, r is the lattice analogue of the continuum 
parameter r. Within our range of parameters and with our choice for t/£ g we found that 
it is enough to set n = 30 in order to eliminate a visible dependence on the two factors 
/ and n individually, making only fn relevant. This ensures that when extrapolating to 
the continuum limit by increasing £g we eliminate sensitivity on the discretization of the 
smearing process. 

The set of matrices g(x) is generated by Monte Carlo for the action Sl- We set £g 
and obtain the corresponding b. We use a combination of Metropolis and Over-Relaxation 
at each site x, where we explore the full SU (N) group. Thus, the operation count goes as 
iV 3 for generating configurations. We found that 200-250 passes suffice to thermalize the 
system starting from g(x) = 1. When going from a configuration with some £g to another 
with £g increased by 1, 50 passes are enough for equilibration and approximate statistical 
independence. 



One major difference between gauge theories [12, 13] and the PCM in the context of 
the large N limit is that in the PCM there is no special suppression of finite volume effects. 
From previous numerical work on the PCM we know that keeping L/^q > 7 reduces finite 
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volume effects to below one percent. As already mentioned, we carry out simulations in 
the range 11 < £g < 20. Therefore we chose L = 150, a relatively large box-size. For 
smearing we chose 

f 2 

t = ^- (5.13) 
300 V ' 

This is a small amount of smearing, but some smearing is necessary in the continuum limit 

because we need to avoid the short distance singularity one would otherwise expect. A 

singularity is incompatible with the W matrices being unitary, and hence bounded, in the 

continuum limit, 

For each g(x) we generate we compute matrices 

W lx {d)=g n {x)gl(x + d- l i) (5.14) 

for positive and negative directions /x and all distinct locations. For a fixed our objective 
is to find d c , the distance at which the eigenvalue distribution of the W's is critical. For 
d < d c the eigenvalue distribution is zero at -1, while for d > d c it is strictly positive 
everywhere on the unit circle. d c is fractional, found by interpolation from the discrete d 
data. The main objective is to show that the quantity 4^ has a nontrivial finite continuum 
limit and to compute that limit. This would establish the transition as a feature of the 
continuum PCM. 

6. The main observable [0. 

As mentioned already, our main observable is related to the average characteristic polyno- 
mial Qn(z,X). A closely related quantity, F(y), is defined by 

F(y, d) = (det( e y/ 2 + e-y' 2 W)) = C + C x y + C 2 y 2 + C 3 y 3 + C 4 y A + . . . (6.1) 

We have 

F(y) = det(W)F*(-y) = F*{-y) = F(-y) (6.2) 

The main reason to introduce F is that it is even in y on account of CP invariance. This 
sets d = C 3 = ... = 0. 

To detect the location of the transition at infinite N we introduce a derived quantity, 
Q(d), which is sensitive to the behavior of F at y = 0. y = corresponds to z = —1 of 
Qn(z, A), which is the point at which the spectrum of the unitary matrix just closes a gap. 



OW/) (6.3) 



£1 is insensitive to overall rescalings of F and y independently. In the DO universality class 
f2 will jump from being constant at | to being constant at | at infinite N. This jump will 
occur at d = d c . The jump in O(ci) can be smoothed out if the critical regime d ~ d c is 
dilated as follows: 

d-d c = ^= (6.4) 
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If we take the correlated N — » do, d — > d c limit, keeping r/ constant and of order unity, we 
obtain a smooth function of r]. At r/ = 0, the value of Q is a universal number: 

r 4 (!■) 

n c = -^Af « 0.365 (6.5) 
We now define a finite N approximation to d c , which we denote d c (N), by 

tt(d c (N),N) = O c (6.6) 
d c without an N dependence is the limit d c (N = oo). One expects: 

d c (N) = d c + ^ + (j±^ (6.7) 

In our simulations we obtained numerical values for £l(d,N) for seven values of d 
centered at 0.65£g and spaced by 1. Using cubic spline interpolation we found d c (N) and 
also an approximation to the derivative of 0(d, N) with respect to d at the point d = d c (N) 
where Q = Q c . We found that this method produces numbers of a relative accuracy of 
fractions of a percent for d c (N) and a few percent for the derivative. Every input set of 7 
fi-values for fixed £g and d was obtained by averaging the coefficients Ck over all possible g- 
pairs and subsequently averaging the Q-s obtained from these average values of Ck for fixed 
configurations of g over about 50 different and statistically independent ^-configurations. 
The errors in the f2-s were obtained by jackknife with single elimination. 

Henceforth, we shall redefine Q to be Vl — Q c . Hence d c (N) will be associated with the 
vanishing of the new Q. 



7. Results. 

For each £q in the range of 11 to 20 we end up with a d c (N) and a slope of Q at £1 = 0. To 
check for a continuum limit we look at the ratios For each value of N we extrapolate 

to the continuum by fitting the data for dc ^ N ' to a second order polynomial in -pr. This 
extrapolation is seen to be quite smooth and the second order polynomial provides a good 
fit. The slopes on the other hand have too low an accuracy to permit an extrapolation to the 
continuum limit; within their errors of a few percent they turn out to be £g independent. 
We took for the slopes the value on the finest lattice we have, at £g = 20 and carried out 
no continuum extrapolation. 

Finally, we are ready to look at the iV dependence. For d c (N) we expect a convergence 
to an infinite iV linear in jj. The slope of d(N) with respect to at ft = should go to 
zero as Sample figures are below. 

Figure [l] shows how the data for is interpolated to get the critical size for a particular 
lattice spacing (the example is for £g = 20) at various values of N. Note that Q gets steeper 
as N increases. The data points markers are larger than the statistical jackknife errors. 

Figure ^ shows the extrapolation to the continuum for iV = 30. The discrete value are 
plotted as squares to stand out, but the sizes of the squares do not indicate the errors. The 
error on the extrapolated value was computed by generating artificial noise on the input 
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Figure 1: Plot of the subtracted as a function of £ for N — 15, 20, 25, 30, 40 and for the finest 
lattice, with = 20. The data point sizes are much larger than the statistical error bars which 
might be visible after significant magnification. 



data with Gaussian errors of size determined by the jackknife estimates of the configuration 
averages. 40 such fake sets were next interpolated by cubic spline and then extrapolated 
to continuum to generate the error on the extrapolated value shown in the figure. 
In Table |T| we summarize our findings. The definition of the entry "slope" is: 

-d(d c (N)/&) 



dtt 



(7.1) 



n=o 



The final objective is to extrapolate to infinite N. 
For the critical size this is shown in Figure We see 
that a simple linear extrapolation in 
obtain in the infinite N limit 



jr works well and 



d c 



continuum; r- 



300 



0.885(3) 



(7.2) 



N 


V ^ G J continuum. 


slope 


15 


0.650(3) 


-4.00 


20 


0.707(2) 


-3.58 


25 


0.741(2) 


-3.22 


30 


0.768(2) 


-2.99 


40 


0.797(2) 


-2.60 



Table 1: Input to the extrapola- 
tion of N to oo 



Even if we accept that the d c (N) that we defined has 
a limit as described, we have presented so far little hard 
evidence that the transition is of the DO type. To see 

this, we look at the slope and check whether numerically it is plausible that it exhibits a 
dependence on y/N as the DO universality class would predict. In Figure ||| we show a plot 
of log |siope| versus log(iV). One sees a slope that is different from the expected -0.5, but 
we have not included any further subleading corrections in -h. There is little doubt that 
a slope of -1 is excluded by the data and the data is consistent with a slope of -0.5 in the 
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Figure 2: Extrapolation to continuum of for N=30. 
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Figure 3: Extrapolation of Icontinuum to infinite N. 



asymptotic large N regime. If we drop the lowest N value from the fit the slope increases 
somewhat, but does not reach 0.5. The slopes for finite N had too little accuracy to allow 
a continuum extrapolation and therefore we used only the slopes at £g = 20, our largest 
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Slope dependence on N 
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Log(N) 

Figure 4: Testing for the critical exponent 1/2 in the slope variable. 

correlation length. There is no reasonable way to get a real error on our slope result, so 
no errors were quoted. Still, our result supports the expectation that the critical exponent 
associated with d — d c is 1/2. 

Nothing we have done so far can be a substitute to a simple direct comparison of the 
DO eigenvalue distribution to the data. In Figure [5] and Figure ^ we show for N=40 plots 
of the eigenvalue density we obtained at £g = 20 for d = 10, 21 and d = 17 respectively. 
Superposed on the data in both figures are the DO distributions obtained by fitting their 
single free parameter k. k gives a measure for how far we are from the critical point which 
corresponds to k c = 2. Durhuus and Olesen obtained an expression for the distribution of 
the angle 9 j§. The eigenvalues of W are e j , with J2j ®j = 0{mod)2iT and p(9)d8 is the 
probability that any 9j satisfy \6j — 9\ < d9/2. p(9) is defined from a complex function 
h(k,9) which depends parametrically on k: 

p(6) = -^Imh(k,0). (7.3) 

The complex valued function h(k, 6) is determined by a nonlinear equation. One 
first introduces two 9 dependent real variables, x and y, and the expressions f(y) = 

isfey g(y)=ycoth( y ). 

x = -f(y)sm(x + 9), x + 9 = Re h(k, 9), y = lmh(k,9) (7.4) 

These equations change variables from h, 9 to x, y, 9. The dependence on k comes in 
through the main equation 

d(y) - k = f(y) cos(x + 9) (7.5) 
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k determines the support of p(9), which just reaches 9 = ±ir at k = k c . 

J 7T if k > 2 

In Figure || we see two extreme cases, one corresponding to a gapped eigenvalue dis- 
tribution for a distance of the order d = 0.5£g an d the other for an eigenvalue distribution 
covering the entire unit circle, for a distance d = 1.05£g- The solid lines are the exact 
DO distribution, at a A; that was determined by the best least squares fit to the empirical 
histogram. Taking data for a larger set of distances then usual, we can ask which distance 
provides a k closest to critical. It turns out that this happens at d = 17, where k = 1.98. 
Thus, this simple method would have given us an estimate for the infinite N critical ratio 
|^ of about 0.85. Our result for our definition for a finite N critical distance d c (N) for 
N = 40 gave ^ ~ 0.80 at £g = 20 while the continuum infinite N number was 0.885. This 
is all consistent. Figure ^ shows this close to critical distribution. 

One may speculate that the DO distribution matches the continuum infinite iV distri- 
bution exactly in the critical case, but this is very likely not correct. Universality holds 
only for 9 close to ±7r. If one looks at the integrated square distance between the DO dis- 
tribution and the empirical one as a function of N there is no evidence that this distance 
extrapolates to zero. 
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N=40, ^ G =20 
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Figure 6: An example of an almost critical eigenvalue distribution (d ~ 0.85£g)- 



8. Summary. 

There is little doubt that with the introduction of the smearing parameter r the PCM 
undergoes an infinite TV phase transition of the DO type, in a manner analogous to smeared 
Wilson loops in two, three and four dimensions. 

The PCM also offers the hope to establish this transition by analytic means. If this 
could be done, the role of the parameter r could be further elucidated. We look at it as 
an extra renormalization of the two point matrix needed to reconcile its short distance 
behavior with it being a fluctuating unitary (and hence bounded) matrix. 

At infinite A" the PCM model has a a trivial S-matrix for its set of bound states. 
Nevertheless, it is not a free field theory in terms of the variables g(x) as evidenced by 
the short distance singularity of the two point function. The introduction of r eliminates 
this divergence but since r can be made very small, the short distance behavior of the two 
point function can be seen for y/r « d « £g- The DO transition occurs at d ~ £g an d 
separates the two regimes, one of noninteracting bound states and the other of Gaussian 
fluctuations of the fundamental variable g(x). 
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